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Abstract 

We consider decay of metastable topological configurations such as strings and do- 
main walls. The transition from a state with higher energy density to a state with 
lower one proceeds through quantum tunneling or through thermally catalyzed quan- 
tum tunneling (at sufficiently small temperatures). The transition rate is calculated at 

^ ■ zero temperature including the preexponential factor and also at a finite low tempera- 

s' ■ 

ture. The thermal catalysis factor is closely related to the probability (effective length) 
of destruction of the string (the domain wall) in collisions of the Goldstone bosons, 
corresponding to transverse waves on the string (wall). We derive a general formula 
which allows to find the probability (effective length) of a string (wall) breakup by a 
collision of arbitrary number of the bosons. We find that the destruction of a string 
only takes place in collisions of even number of the bosons, while the destruction of 
the wall can occur in a collision of any number of particles. We explicitly calculate 
the energy dependence of such processes in two-particle collisions for arbitrary relation 
between the energy and the largest infrared scale (the size of a critical gap). 
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1 Introduction 

Topological configurations of fields with the geometry of a string or a domain wall arise in 
various models either as solutions of classical field equations or as nonperturbative effective 
configurations of essentially quantum fields. String-like configurations are present in poly- 
mers, superconductors, in theoretical models of non-Abelian dynamics, as well as in models 
of QCD confinement, and as topological defects in models with spontaneous breaking of 
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gauge or global symmetries. Metastable domain wall solutions arise in models with spon- 
taneously broken approximate symmetry. The existence of such solutions was shown from 
different points of view, for example, in Refs. [1, 2, 3]. 

These configurations are classically stable, since there is no classical trajectory which 
takes from one topological solution into another with different topological number. But it 
may be the case that a classically forbidden trajectory for such a deformation exists so that 
the considered configurations are unstable quasiclassically and decay due to e.g. quantum 
tunneling effects. As an illustration of such decay process one can consider the example of 
a metastable domain wall investigated in Refs. [4, 5] which arises in a model with a scalar 
field potential shown in Fig.l. The domain wall solution corresponds to the interpolation 
between the same vacuum state at two spatial infinities, e.g. at z — — oo and z = +oo with 
the field winding around the 'peg' in the potential. Such configuration is classically stable 




Figure 1: Potential 



in the sense that the solution with certain winding number can not evolve classically to a 
solution with different topological number. However, such a transition can proceed either 
due to a thermal fluctuation, where the path, corresponding to the solution, is lifted over 
the barrier or due to quantum tunneling (the motion in a classically forbidden region). 

The topological configurations can be metastable with respect to either a complete break- 
ing, or a transition to an object of lower tension emerging instead of the initial one. The 
former situation is relevant e.g. for a break up of a QCD string with formation of a quark 
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- antiquark pair, or a formation of a monopole - antimonopole pair [6], and also in a whole 
class of theories with spontaneous symmetry breaking [7]. The latter situation involving a 
phase transition between states of a string with different tension is found e.g. in Abelian 
Higgs models embedded in non- Abelian theories [8, 9]. 

A domain wall decay is analogous to the well-known spontaneous false vacuum decay 
process [10, 11, 12] in 2 + 1 dimensions. Indeed, if a hole of area A c is created in a wall with 
tension e, the gain in the energy is eA c - The barrier that inhibits the process is created by 
the energy aV, with V being the perimeter of the hole and a being a tension associated with 
the interface. Thus, the 'area' energy gain exceeds the barrier energy only starting from a 
critical size of the hole created, i.e. starting with a round hole of radius R w = 1o je. Once a 
critical bubble of the lower phase has nucleated due to tunneling, it expands, converting the 
domain wall. Therefore the probability of the transition is given by the rate of nucleation of 
the critical holes in the domain wall. 

The same argument also applies to the string decay which is similar to a metastable 
vacuum decay process in 1 + 1 dimensions. The critical size of the gap in the initial string 
is then i c = 2R S = 2fi/e. Moreover due to the geometry of the problem such decay also 
bears a great similarity to the well known Schwinger process of production of pair of charged 
particles by an external electromagnetic field [13]. 

Despite the similarity with the false vacuum decay the breaking of the strings and walls 
involves an essential difference associated with the transverse waves propagating on the 
latter objects, corresponding to massless bosons, which in fact are the Goldstone bosons of 
the spontaneously broken translational invariance. The effect of the soft modes of the field of 
these bosons enters the probability of breaking the topological defects at the preexponential 
level [5, 9]. 

The process of a metastable string transition from a state with tension E\ to a state with 
tension e 2 was considered in [6, 7, 8]. Using the analogy with the false vacuum decay the 
decay rate was found with exponential accuracy in terms of the rate of nucleation of critical 
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gaps per unit length of the string, 7 S = dT/di 1 , 

ls = Ce^(--^—\ , (1) 

where \x is the mass associated with the interface between the two phases of the string. The 
decay rate of an axion domain wall was considered at the level of the semiclassical exponent 
in Ref. [4] by adapting the expression for metastable vacuum decay in 2 + 1 dimensions. 

The preexponential factors for the spontaneous decay of string and walls were calculated 
in [5, 9]. For a string the rate of the transition between the states with tension e 1 and e 2 
was found to be 



Is ~- 



F 



2tt 

with the dimensionless factor F given by 




(2) 



F m = x ^±^T( £ -l±^ + l) ( e J~Il)^ exp ( £ -l±^) (2, £ -l±^X 1 ' 2 , (3) 
\Ei/ V 2e i / \e 1 + e 2 J \e 1 -e 2 J\ e 1 — e 2 J 

while the decay rate of a domain wall is formulated in terms of the rate ^ w of nucleation of 
critical holes per unit area, 7 W = dT/dA, 

where hr and an are renormalized mass and tension parameters associated with the interface 
of the strings and domain walls respectively, and C is a constant that does not depend on e 
but does depend on other dimensional parameters in the underlying field theory. 

The thermal effects in the decay of metastable strings [14] further expose the difference 
from the false vacuum decay and the Schwinger process. Namely, as long as the temperature 
T is lower than the inverse of the critical length, T < 1/4 the thermal effects in the latter 
processes are exponentially suppressed as exp(— m/T) with m being the lowest scale for 
particle masses in the theory, and these effects are very small due to the strongly suppressed 
presence of massive particles in the thermal equilibrium [15]. In the case of a string, however, 



x It should be mentioned that in d — 2 there is no transverse motion for the string, and the preexponential 
factor can be also copied from the known result of spontaneous metastable vacuum decay 

Gd = 2 ~ ■ 
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the transverse waves on the string are massless so that their excitation has no suppression 
by the mass at arbitrarily low temperature. The thermal excitations of these waves create 
fluctuations in the distribution of the energy in the string which catalyze the nucleation of 
the critical gap. Clearly, at T <C l/£ c the typical wavelength of the thermal waves ~ 1/T 
is large in comparison with the critical length £ c , and the thermal effect in the rate is quite 
small, although not exponentially small. The leading low temperature correction in the 
nucleation rate is given by the thermal catalysis factor[14] 



while as T approaches 1/4, the catalysis factor develops a singularity at £ C T = 1. At still 
higher temperatures the considered string transition behaves, in a sense, similarly to the 
false vacuum decay [15], namely the regime of the transition changes to a different tunneling 
trajectory, so that the temperature dependence appears in the semiclassical exponential 
factor in Eq.(2), rather than in the preexponential term. 

The effect of the thermally excited waves in the decay rate can in fact be considered [16] 
as an additional contribution to the probability of the decay due to collisions of the Goldstone 
bosons that are present in the thermal bath. As it turns out, it is possible to identify in 
the thermal catalysis factor the contribution of individual such collisions and thus calculate 
the probability of the destruction of metastable string by colliding particles. In particular, 
the first temperature dependent term in the expansion (5) is entirely due to the process 
of the critical gap creation in a collision of two particles in the limit, where their center 
of mass energy E = y/s is much smaller than l/£ c . In order to separate the terms in the 
decay probability originating from collisions of different number of the bosons, the standard 
thermal approach to calculating the decay rate at finite temperature [14] is modified [16] by 
formally introducing a negative chemical potential for the bosons. This allows to find the 
dependence of the probability of string destruction in an n-boson collision for an arbitrary 
relation between £ c and the energy of the bosons. One of the results of such a consideration 
is that the string destruction takes place only in collisions of even number of bosons and 
is absent at odd n. The explicit expression for the (dimensionless) probability W 2 of the 
critical gap creation in two-boson collisions at arbitrary values of the parameter E £ c has 
especially simple form for the decay of the string into 'nothing' i.e. at e 2 = 0: 






(6) 
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with Is(x) being the standard notation for the modified Bessel function of the third order. 
This expression is applicable as long as the energy E is small in comparison with the string 
energy scale: E <C y 7 ?, which condition still allows for the parameter E £ c to be large. At 
E £ c ^> 1 the energy dependent factor in Eq.(6) has the exponential behavior exp(E£ c ) which 
matches the semiclassical approximation for the energy-dependent factor [17], corresponding 
to tunneling at the energy E rather than at zero energy. In this sense the discussed en- 
ergy dependence of the collision-induced probability of the decay describes the onset of the 
semiclassical behavior. 

One can also notice that the thermal factor IC S , is singular [14] at £ C T = 1. However the 
two-boson production rate does not exhibit any singularity at any value of the parameter 
E£ c , and a similar smooth energy behavior is also true for individual n-boson processes. 
Thus the 'explosion' of the thermal rate at £ C T = 1 is a result of infinite number of processes 
becoming important at this point, rather than due to a finite set of processes with a limited 
range of n developing large probability at the energy per particle of the order of \ji c . 

The purpose of this paper is to generalize to the case of a metastable domain wall the ap- 
proach used for analyzing the thermal and collision induced effects in the decay of metastable 
strings. We find the general expression for the thermal catalysis factor K, w . In particular, 
the leading term in this factor at low temperature (in d dimensions) is given by 

K w = 1 + 12 (d - 3) C 2 (5) (R W T) W + — (7) 

where R w is the radius of the critical hole in the metastable wall, and ((x) is the standard 
notation for the Riemann £ function, £(5) ~ 1.037. We also aim at providing here a sys- 
tematic and self-contained presentation of the calculations leading to our final results. To 
this end we recapitulate in detail in Section 2 the technically simpler calculations for the 
processes with metastable strings, and then in Section 3 expand the method to the case of 
a metastable domain wall. Such 'semi- review' format of the paper also allows us to illus- 
trate the similarity and the difference between the string and the domain wall transitions. 
Although the calculations in both cases are very much similar, the results differ in some 
essential details. In particular the destruction of a metastable wall is induced by a collision 
of any number of Goldstone bosons, whereas a string can be destroyed by a collision of only 
even number of the bosons. Another difference relates to the thermal catalysis of the decay: 
for a string the low temperature expansion for the enhancement of the quantum tunneling 
is applicable everywhere where the series is convergent, i.e. up to the temperature T = l/£ c , 
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while for a domain wall the thermal fluctuations start dominating [15] at a lower tempera- 
ture (T = 3/(4£ w ) with i w being the diameter of the critical hole in the wall), at which the 
thermal effects in the quantum tunneling probability are still very small numerically. 

2 String-like configuration 

2.1 Spontaneous decay of a metastable string 

The tunneling trajectory can be described in the Euclidean space by a configuration called the 
'bounce' [11], which is a solution to (Euclidean) classical equations of motion. The general 
expression for the effective Euclidean space action for the string with the two considered 
phases can be written in the familiar Nambu-Goto form: 

S = fiP + e 1 A 1 +e 2 A 2 , (8) 

where A\ and A 2 are the areas of the world sheet for the two phases, and P (the perimeter) 
is the length of the world line for the interface between them. 




Figure 2: The bounce configuration, describing the semiclassical tunneling trajectory. 

The action (8) is an effective low-energy expression in the sense that it only describes the 
'stringy' variables and is applicable as long as the effects of thickness of the string and of 
its internal structure can be neglected. Denoting Mo the mass scale at which such approach 
becomes invalid (e.g. the thickness of the string r ~ 1/M ), one can write the condition 
for the applicability of the effective action (8) in terms of the length scale, i ^> 1/M , and 
the momentum scale k <C Mq. Assuming that the initial very long string with tension e± is 
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located along the x axis, one can readily find that the action (8) has a nontrivial stationary 
configuration, the bounce, namely, that of a disk in the (t,x) plane occupied by the phase 2, 
as shown in Figure 2, with the radius 

R s = , (9) 

ei - e 2 

which is the radius (one half of the length £ c ) of the critical gap. The difference between 
the action (8) on this configuration and on the trivial one is exactly the expression for the 
exponential power in Eq.(l), and the condition for applicability of the effective action (8) 
requires 

M R s = ^^»l. (10) 
£i - £2 

Generally one also has /i>M , and for the strings in weakly coupled theories \i ^> M , so 
that the power in the exponent in Eq.(l) is large, which justifies a semiclassical treatment. 

The probability of the transition is determined [11, 18, 19] by (the imaginary part of) the 
ratio of the path integrals Z\ 2 and Z\ calculated with the action (8) around respectively the 
bounce configuration and around the initial flat string: 

?. = ilmf. (ID 

It can also be reminded that, as explained in great detail in Ref. [18], that the imaginary 
part of Z\2 arises from one negative mode at the bounce configuration, and that due to two 
translational zero modes the numerator in Eq.(ll) is proportional to the total space time 
area A in the (t, x) plane occupied by the string, so that the finite quantity is the transition 
probability per unit time (the rate) and per unit length of the string. 

In order to evaluate the relevant path integrals with pre-exponential accuracy we use 
the cylindrical coordinates, with r and being the polar variables in the (t, x) plane (of the 
bounce), and z being the transverse coordinate. We consider only one transverse coordi- 
nate, since the effect of each of the extra dimensions factorizes, so that the corresponding 
generalization is straightforward. We further assume, for definiteness, that the space-time 
boundary in the (t, x) plane is a circle of large radius L, where the boundary condition for 
the string is z{r — L) — 0. The small deviations of the string configuration from the bounce, 
illustrated in Fig.3, can be parametrized by the radial (/) and transverse (£) shifts of the 
boundary between the string phases: 

r{9) = R S + f(9) , z{9) = £(0) , (12) 
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and by the variations of the surfaces of the two string phases: zi(r,0) and 2:2(7*, 9), where, 
naturally, 

z 1 (R s ,9) = z 2 (R s ,9)=£(9) . (13) 



€ 2 




Figure 3: Fluctuations around the bounce configuration. 



In terms of these variables the action (8) can be written in the quadratic approximation 
in the deviations from the bounce as 

S 12 = e 1 nL> + ^ + £ -±^ Ue (e 2 + / 2 -/ 2 ) + 
£1 - £2 2 7^ 1 

t/> rf 4 ;2+ ^Hf r ^ 2+ ^)- (i4) 

where the primed and dotted symbols stand for the derivatives with respect to r and 9 
correspondingly. 

Finally, the action around a flat initial string configuration in the quadratic approxima- 
tion takes the form 

S 1 = £ 1 vrL 2 + - [ L 'rdrdO iz' 2 + , (15) 
2 Jo \ r 2 ) 

with z(r, 9) parametrizing small deviations of the string in the transverse direction. 
2.1.1 Separating variables in the path integrals 

One can readily see that in the quadratic part of the action (14) the 'longitudinal' variation of 
the bounce boundary in the (t, x) plane, described by the function f{9) completely decouples 
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from the rest of the variables. This implies that the path integral over / can be considered 
independently of the integration over other variables and that it enters as a factor in Z 12 . On 
the other hand it is this integral that provides the imaginary part to the partition function, 
and it is also proportional to the total space-time area A. Moreover, this path integral is 
identical to the one entering the problem of false vacuum decay in (1+1) dimensions and we 
can directly apply the result of that calculation [20]: 



£1 ~ £2 



(16) 



2tt 

The expression for the transition rate thus can be written in the form 

-77 = ex P , ( 17 ) 

at 2n \ e 1 — e 2 J Z\ 

with the path integral Z 12 running only over the transverse variables £, z\ and z 2 , 

Z12 = J ViVz x Vz 2 exp (~S 12 ) (18) 
and involving only the quadratic in these variables part of the action (14), 

5l2 = £ip/ dsf + |£ r ^( z ; 2 + ^ + |/;V^( 42 + ^ . (19) 

In the same quadratic approximation the flat string partition function Z 1 is given by 

Z 1 = Jvz exp (-Si) (20) 

with Si given by Eq.(15) and the integral running over all the functions vanishing at the 
space-time boundary: z(L, 6) = 0. 

At this point there still is a coupling in the path integral (18) between the bulk variables 
Zi, z 2 and the boundary variable £ arising from the boundary conditions (13). This however 
is a simple issue which is resolved by a straightforward shift of the integration variables Z\ 
and z 2 . Namely, we write 

zi(r, 9) = z lc (r, 9) + z lq (r, 9) , z 2 (r, 9) = z 2c (r, 9) + z 2q (r, 9) , (21) 

where z\ q and z 2q are the new integration variables in Z\ 2 and these functions satisfy zero 
boundary conditions, 

z lq (R s , 9) = z 2q (R s , 9) = z lq (L, 9) = z 2q (L, 9) = , (22) 
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while z\ c and z 2c are fixed (for a fixed £(0)) functions satisfying the boundary conditions 

z lc (R s , 6) = z 2c (R s , 6) = S(6) , z lc (L, 6) = , (23) 

which are also harmonic, i.e. satisfying the Laplace equation Az = 0. 

After the specified shift of the variables we arrive at the expression for the action (19) in 
which the bulk and the boundary degrees of freedom are fully separated: 



Sl2 — 



J d0i 2 + R s J d0 {^d r z 2c \ r=Rs - y&Zlc \ r=R ^j £ 



- y / <Prz lq Az lq - y J d 2 rz 2q Az 2q . (24) 

Clearly, the boundary terms in the first line here, arising from the integration by parts, 
depend only on the transverse shift of the boundary £(#). The partition function Z u can 
thus be written as a product of the 'boundary' and the 'bulk' terms: 

Z 12 = ^12(boundary) ^12(bulk) , (25) 

with the 2i2(buik) being given by the path integral over the bulk variables z\ q and z 2q only : 
-Zi2(buik) = J Vz lq Vz q exp ^y J d 2 r z lq Az lq + y j d 2 r z 2q Az 2q ^j , (26) 
and the boundary term 

-2l2(boundary) = j V£ exp - — j d0 £ 2 - R s J d6 {^d r Z 2c \ r=Rs - ^8 r Z lc \ r=R )j £ 

(27) 

involving integration over only the boundary values. 

The subsequent calculation of the ratio of the partition functions in Eq.(17) can in fact 
be reduced to a calculation of i?i2(boundary) only. In order to achieve this reduction one should 
organize the partition function Z\ in the denominator of Eq.(17) in a form similar to Eq.(25) 
as follows. Although the flat string configuration makes no reference to a circle of the radius 
R, the partition function Z\ can be calculated by first fixing the transverse variable z at 
r = R s : z(R s ,0) = £(6) and separating the integration over the bulk variables. Then the 
flat string partition function factorizes in the form similar to Eq.(25): 

%1 = %1 (boundary) ^l(bulk) j 

(28) 

with -Zi(buik) given a similar path integral as in Eq.(26), 

^i(buik) = f Vz lq Vz q exp ^y J d 2 r z lq Az lq + y j d 2 r z 2q Az 2q ^j , (29) 
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where, as in Eq.(26), z\ q and z 2q are respectively the outer (i.e. at r > R s ) and the inner (r < 
R s ) transverse fluctuations with zero boundary conditions. The difference in the coefficient in 
the expressions (26) and (29) for the contribution of the inner part, e 2 vs. e±, is not essential, 
since the overall coefficient of the quadratic part of the action is absorbed in the measure of 
integration, as can be seen by rescaling to the corresponding canonically normalized variables 
= y/ez. 

One therefore finds that -Zi(buik) — -2^i2(buik), and the ratio of the partition functions in 
Eq.(17) is in fact determined by the ratio of the boundary terms. 

2.1.2 Regularization 

The boundary factor 2i(boundary) for the flat string is somewhat different from the one given 
by Eq.(27) and reads as 



Z\ (boundary) = J T>£ exp -i? s | J d6 (d r z 2c \ r=Rs - d r z lc \ r=Rs ) € 



(30) 



where the functions z\ c and z 2c are defined in the same way as in Eq.(27). 

The latter functions can be readily found by expanding the boundary function £(#) in 
angular harmonics: 

£(6) = -L + — ]T [a„ cos(n^) + b n sin(n0)] (31) 

V 27T y/K n=l 

with a n and b n being the amplitudes of the harmonics. Then at n ^ the harmonics of the 
discussed functions are found as 

1 R n 1 v n 

Q) = —i= K cos{n6) + b n sin(n0)] -± , z£ ] (r, 0) = ~f= [an cos{n6) + b n sin(n0)] — , 

y/TT r n y/TT tig 

(32) 

and for n=0 these are given by 

Substituting these expressions for the harmonics in the equations (27) and (30) and perform- 
ing the Gaussian integration over the amplitudes a n and b n we find the boundary factors in 
the form 

^Kb—rt = ft ^ In - n (fl _ £2) , i2 + (£l + £2) „ (34) 
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and 



-Zi(b. 



oundary) 



N \ Ex In 



n 



(35) 



with Af being a normalization factor. 

Clearly, each of the formal expressions (34) and (35) contains a divergent product, and 
their ratio is also ill defined, so that our calculation requires a regularization procedure that 
would cut off the contribution of harmonics with large n. A regularization of high harmonics 
is also required on general grounds. Indeed, as previously mentioned, our consideration 
using the effective string action (8) is only valid for smooth deformations of the string, i.e. 
as long as the relevant momenta are smaller than the mass scale M for excitation of the 
internal degrees of freedom within the thickness of the string. For an n-th harmonic the 
relevant momentum is k ~ n/R s so that the applicability of the effective low energy action 
requires a cutoff at n <C M R S . In order to perform such regularization we use the standard 
Pauli-Villars method and introduce a regulator field Z with negative norm and the action 
corresponding to the quadratic part of the Nambu-Goto expression (8) for small z: 



1R 



2 J A x 



(d^Z) 2 + M 2 Z 2 ] + — / d 2 r Ud^Z) 2 + M 2 Z' 
J 2 Ja 2 l 



(36) 



with M being the regulator mass, which physically should be understood as satisfying the 
condition M M and still being much larger than the relevant scale in the discussed 
problem, in particular MR S ^> 1. 

The regularized expression for the ratio of the boundary terms in Zyi and Z\ thus takes 
the form 



Z 



12(boundary) 



K 



^12(boundary) 
12(boundary) 



-2l (boundary) 



Z\ (boundary) 

where we introduced the notation 1Z for the regularized ratio, and the regulator partition 



MR) 
1 (boundary) . 



(37) 



functions Z 



(R) 

12(boundary) 



and Z 



(R) 

1 (boundary) 



are determined by the same expressions as in Eqs.(27) 



and (30) with the 'outer' and 'inner' functions Z\ c and z 2c being replaced by their regulator 
counterparts Z\ c and Z 2c which still satisfy the boundary conditions similar to (23): 



Zic(R s ,0) — Z 2c (R s ,9) — £r(6) , 



(38) 



but are the solutions of the Helmholtz rather than Laplace equation: (A — M 2 )Z = 0. 

The solutions of the latter equation fall off exponentially at the scale determined by M, 
and for our purposes only the behavior near the circle r = R s is needed. For this reason we 



13 



write the equation for the radial part of the n-th angular harmonic Z n (r), 

,2 



Zn+ 1 -Z' n -^Z n -M^Z n 



, (39) 

and parametrize the radial coordinate as r = R s + x, and treat the parameter (x/R s ) as 
small, since the scale for the variation of the solution is x ~ 1/^M 2 + n 2 / R 2 S . This approach 
yields an expansion of the regulator action associated with the boundary at r = R in powers 
of 1/ \J (MR S ) 2 + n 2 . With the accuracy required in the present calculation, the (normalized 
to one at r = R s ) solution to Eq.(39) is found in the first order of expansion in (x/R s ) as 



Z n (R s + x) 



1 - 



(MR S 



x 



2 (MR S ) 2 + n 2 R s 



exp -^J(MR S ) 2 + n 2 



\x\ 
R.« 



(40) 



Using the form of the solutions for the harmonics of the regulator field given by Eq.(40) 
and the expressions (34) and (35), one can write the regularized ratio of the boundary 
partition functions (37) as 



K 



lei+e 2 
2e 1 



n 

n=l 



n 2 + bJ(MR s ) 2 + n 2 



n 



n=l 



1 + 



n 2 + bn 

(MR,) 2 



TT r l 

„=1 ^(MR S ) 2 + 



X 



2 [(MR S ) 2 + n 2 } [n 2 + b yj (MR S ) 2 + n 2 
where we introduced the notation 



(41) 



(42) 



and the last product factor in Eq.(41) arises from the first term of expansion in (x/R) in the 
pre-exponential factor in Eq.(40) 

2.1.3 Calculating the products 

Each of the products in Eq.(41) is finite at a finite M and can be calculated separately. We 
start with the most straightforward one, which is directly given by an Euler's formula: 



n 



n 



tiMR, 



n=i ^J(MR S ) 2 + n 2 V sinh( 7 rM J R s ) 
where the last transition corresponds to the limit MR S ^> 1. 



2nMR s exp 



7iMR< 



(43) 
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The other two factors in Eq.(41), the first and the last, generally depend on the relation 
between the parameters b and MR S , or equivalently between {e\ + £2) and fiM. We find 
however that the latter product is equal to one at MR S ^> 1 independently of b. In particular 
in the nontrivial case of b <C MR S we find 



inn 

71=1 



1 + 



(MR s f 



2 [(MR s ) 2 + n 2 } n 2 + bJ(MR s ) 2 + n 



[ 2 J no 



{MR s ) 2 + n 2 



n 2 + b ^J(MR S ) 2 + n 



l -11 



, (44) 



MRs 



where the lower limit in the integral is any number uq that is finite in the limit MR S — > 00. 

The dependence of the first product factor in Eq.(41) on the ratio {e\ + 82) /{fiM) = 
b/(MR s ) is essential and we consider two limiting cases when this ratio is much bigger than 
one and when it is very small. In the former case, i.e. for b ^> MR S , the first product in 
Eq.(41) becomes reciprocal of the second, and one finds 



()»Mfi s »l 



1 . 



(45) 



(Clearly one can also safely make the replacement (e\ + £2)/(2ei) — > 1 at b 1.) 

The behavior of 1Z in the opposite limit, i.e. at b <C MR S , turns out to be significantly 
more interesting. Using the Euler-Maclaurin summation formula for the logarithm of the 
first product in Eq.(41) we find in the limit MR S ^> 1 and MR S ^> b: 



°° n 2 + bJ{MR s ) 2 + n 2 T(b+1) 
11 ~> — = o_ rrrrrr ex P 



, n 2 + bn 2n^bMR, 

Being combined with the expression (44) this yields the formula 



7T JbMR s - b \n(MR s ) - b (1 - In 2) 



• (46) 



n = 



le 1 +e 2 V(b + i) 



2£i V2V6 



exp 



-- MR S + 7T ^jbMR s - b \n(MR s ) - b (1 - In 2) 



(47) 



which contains an essential dependence on the regulator mass parameter M. We will show 
however that all such dependence in the phase transition rate can be absorbed in renormal- 
ization of the parameter jj, in the leading semiclassical term. 



2.1.4 Renormalization of /x 

The parameter /i is defined in the action (8) as the coefficient in front of the length of the 
boundary between the world sheets for two phases of the string. Generally this parameter 
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gets renormalized by the quantum corrections, and in order to find such renormalization at 
the level of first quantum corrections, one needs to perform the path integration using the 
quadratic part of the action around a configuration, in which the length of the interface is an 
arbitrary parameter. For a practical calculation of this effect we consider a Euclidean space 
configuration, shown in Fig.4, with the string lying flat along the x axis, and the interface 
between the two phases being at x — 0. The length of the world line for the boundary is thus 
the total size T of the world sheet in the time direction. It should be mentioned that such 
configuration with different string tension on each side of the boundary is not stationary for 
the action (8). However it can be 'stabilized' by a source term (external force) depending on 
the coordinate x(t) of the boundary: J J(t)x(t) dt, which term does not affect the quadratic 
in fluctuations part of the action. 



Figure 4: The configuration for the calculation of the renormalization of \i. 

The Gaussian path integral over the transverse coordinates z(x, t) is then to be calculated 
with zero boundary conditions at the edges of the total world sheet. Using the notation £(t) 
for the transverse shift of the boundary, this condition implies £(0) = £(T) = 0, so that the 
function £(£) has the Fourier expansion of the form 



and a similar expansion applies to the regulator boundary function £_r(£). The part of the 
effective action associated with the boundary is determined by the functions z c (x,t) for the 
transverse shift of the string and the corresponding regulator functions Z c (x,t) that satisfy 
the equations 





(48) 



Az c = and (A - M 2 ) Z c = 



(49) 
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and the boundary conditions 



z c (o,t) =t(t) , z c (o, t)=£ R (t) 



(50) 



as well as zero boundary conditions at the edges of the total world sheet. One can readily 
find these functions for each harmonic of the boundary values £ and £ R , using the solutions 
for z r and Z r in each harmonic: 



x,t) = exp ( -\x\ — 



Tin 



sin 



and 



ZW(x,t) 



exp 



\x\ 



Tint 



sin 



Tint 



(51) 



(52) 



In order to separate the boundary effect in the path integral around the considered 
configuration from the bulk effects, we again divide it by the path integral around the 
configuration where the whole world sheet is occupied by the same phase of the string. 
The latter phase can be chosen with either of the tensions, or with an arbitrary tension e. 
Such division results, as previously, in the cancellation of the bulk contributions, and the 
remaining part of the effective action associated with the boundary is written in terms of 
the regularized path integral over the boundary function £ as 

fx R T = fiT - 

/P£exp{-i JdtLe + m (e 2 z / c (x,t)\ x ^ -e 1 z' c (x,t)\ x ^ +0 )]} 

In 7^ F - ^ ; i N + 



/D&expj-i Jdt^ R 2 + Ut) (e 2 Z' c (x,t)\ x ^ -e 1 Z' c (x,t)\ x ^ +0 )^ 



In 



yPj exp{-e J dt^t) z' c (x,t)\ x ^ } 
JVU exp{-e fdtt R (t)Z> c (x,t)\ x ^ } ' 



(53) 



where fj, R — /i + 5/i is the renormalized mass parameter. The correction to /i can thus be 
written in the form 



where 



n 2 + bn t=i ^j{MT/n) 2 + n 2 

£i + e 2 T 



b = 



fl 71 



(54) 



(55) 
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In the limit £1 + 62 <C (i>M one can directly apply the results in Eqs.(43) and (46) for 
evaluation of the expression (54) and write 



1 / ~ ~ ~ MT ~ MT ~ MT ~ 
5fx = --=- 61n6-6-^ + 7ri/6— -61n— -6(l-ln2) 

2 T \ 2 V TT 7T 



M 1 {e 1 + e 2 )M e 1 + e 2 2( £l + £ 2 ) 

= T~2V ^ 2^T ln ^— ' (56) 

where a use is made of the Stirling formula 

to r(5 ±j ) ->5(hi6-i), 



2tt6 

considering that 6 is proportional to large T. 

In the limiting case where £1 + e 2 ^> A*-M the correction 5/j, vanishes, so that the renor- 
malization effect is negligible. 

2.1.5 The result 

We can now assemble all the relevant elements into a formula for the rate of the considered 
transition. Clearly, the path integration over the variables z factorizes for each of the (d — 2) 
transverse dimensions, so that the expression for the decay rate takes the form 

7 s = ^^- 2 exp(-^) , (57) 
2tt \ £i-e 2 J 

where \x is the zeroth order mass parameter. In the case of large string tension, E\ + £ 2 ^> 
//Mo, the 'bare' \x coincides with the renormalized one, and the factor 1Z is equal to one. 
It can be noted that the resulting obvious expression for the rate is also correctly given by 
Eq.(2) as soon as the factor F in Eq.(3) is taken in the limit E\ — e 2 S\ + e 2 : F — > 1, 
which limit, as previously discussed, is mandated in this case. 

In the opposite limit of heavy /x, S\ + e 2 /iM , both the expression (47) depends on 
the regulator mass M and the M-dependent renormalization of /j is essential. Taking into 
account that each of the transverse dimensions contributes additively to 8/j, and expressing 
in Eq.(57) the bare /i through the renormalized one: /x = — 5fi, one readily finds that the 
dependence on the regulator mass M cancels in the transition rate, and one arrives at the 
formula given by Eq.(2) and Eq.(3). 
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The formula (3) is applicable for arbitrary ratio of the string tensions 82/ £\- In particular 
it can be also applied at 62 = 0, in which case the considered transition describes a complete 
breaking of the string. 

It can be also noted that numerically the factor F depends very moderately on the ratio 
of the tensions and changes approximately linearly between -F(O) = e/v^47r = 0.7668 . . . and 
F{1) = 1. We thus conclude that the two-dimensional expression (e\ — 62) /(27r) for the pre- 
exponential factor in the transition rate provides a fairly accurate approximation in higher 
dimensions as well, as long as the exponential factor is expressed in terms of the physical 
renormalized mass /zr. 

There is however an interesting methodical point pertaining to the considered here prob- 
lem. Indeed, as was already mentioned, the difference from the problem of particle creation 
by external electric field is in that the motion of the ends of the string involves in addition 
to the mass /z also an adjacent part of the string. In terms of the calculation of the path 
integrals around the bounce the difference is in that the spectrum of soft modes in the parti- 
cle creation problem (as well as in that of the two-dimensional false vacuum decay) consists 
only of the modes associated with one-dimensional world line of the boundary of the bounce. 
The entire pre-exponential factor can then be found using the effective low energy action for 
these modes[20]. In the considered here string transition there are also low modes in the bulk 
of the world sheet of the string, and there is no parametric separation of their eigenvalues 
from those of the modes associated with fluctuations of the boundary. In the presented cal- 
culation the separation of the boundary and bulk variables is achieved through an 'artificial' 
organization of the normalization partition function for a flat string into boundary and bulk 
factors boundary and ^buik- The bulk contribution then cancels in the ratio of the partition 
functions near the bounce and near a flat string, so that the remaining calculation is reduced 
to considering the integrals over the boundary functions only. One can also readily notice 
that the additional contribution to the action from the boundary terms as e.g. those with 
the functions z\ c and Z2 C in Eq.(27) corresponds to precisely the effect of 'dragging' of the 
string by its end. 

2.2 String transition at non-zero temperature 

The formula in Eq.(ll) corresponds to a calculation of the decay rate as the imaginary part of 
the energy of the initial string. At a finite temperature T the corresponding relevant quantity 
is the imaginary part of the free energy[21], which one can calculate in the Euclidean space 
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by considering periodic in time configurations with the period (3 = T . In other words the 
thermal calculation corresponds to the path integration in the Euclidean space-time having 
the topology of a cylinder. The nucleation rate is then described by the same formula (11) 
with the action and the area being calculated over one period, A = X (3, where X is the 
length of the string along the spatial dimension. 

We consider only sufficiently small temperatures (3 > £ c , at which temperatures we show 
the thermal effects behaving as powers of T, which distinguish the string process from the 
decay of metastable vacuum[15]. We also treat the length X of the metastable string as the 
largest length parameter in the problem, so that /3<X. Under these conditions the bounce 
corresponding to the action (8) is the same as at zero temperature, except that it is placed 
on a long cylinder rather than on a large flat plane (Fig. 5). 



As above we aim at calculating the path integral over the variations of the string around 
the bounce configuration as illustrated in Fig. 6. The coordinates on the cylinder (or on 
the plane where all the points separated by nf3 along Euclidean time are identified) £1X6 t ; 37 . 
with t being the periodic time coordinate, and the coordinate orthogonal to the surface of 
the cylinder is z. The boundary conditions for the configurations over which we integrate 
are 




E 2 



Figure 5: Bounce configuration for nonzero temperature 





(58) 
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In polar coordinates (r, 9) in the (t, x)-plane one finds the variation of the bounce action (8) 
similar to (14) while for the flat string one has the expression analogous to that in Eq.(15). 




E 2 



Figure 6: Variation of the bounce configuration 



As it was shown for the case of zero temperature the partition functions factorize, there- 
fore the expression for a decay rate is given by (17), with the only difference for a nonzero 
temperature case that the boundary conditions for outer solution are 



while the inner solution z 2c is required to be regular inside the disk. 

In order to do the path integrals as before one can expand the boundary function 
in angular harmonics. For the inner solution z 2c to the Laplace equation, i.e. at r < R s one 
finds no difficulty in finding the harmonics matching the boundary function at the interface 
(32). For the outer solution z\ c however there is a difficulty due to the mismatch between 
the symmetry of the boundary and of the periodicity conditions. It is impossible to choose 
the solution to the Laplace equation for outer string bulk variable to be z[™\r, 6) ~ r _n , 
since it is not periodic in time. In this situation in order to have a periodic solution one can 
perform a periodic mapping of the cylinder on the plane and consider the outer solution of 



z lc {t + (3,x) = z lc (t,x), 



(59) 



and 




(60) 
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the Laplace equation as the sum of the solutions produced by a "source" at each period as 
illustrated in Fig. 7. Introducing the complex variable w = t + ix, we construct the solutions 
for the functions Zi c (r, 9) using the harmonic real and imaginary parts of the following basis 
set of periodic functions, satisfying the boundary condition (60) at large \x\. 




Figure 7: Periodic configuration on the plane 

Clearly, these functions are periodic in the Euclidean time by construction, with the 
period (3. Also the functions g k (w) with k > 1 are analytic complex functions, so that 
their real and imaginary parts are harmonic. In the functions go(w) and gi(w) the explicit 
dependence on x, introducing non-analyticity, is linear and is thus also harmonic, so that 
their real and imaginary parts do satisfy the Laplace equation. 

An arbitrary periodic outer solution to the Laplace equation, satisfying the boundary 
condition (60) at large |a;| can be expanded in the series 

oo 

zu(r, 9) = A Reg (w) + J2A k Reg k (w) + B k lmg k (w), (62) 

k=l 
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The disadvantage of this set of solutions (hence of such an expansion) is that it is not 
orthogonal, so that the expression for the action up to quadratic terms is not diagonal in 
this basis. Therefore, the calculation of the integral is not just a calculation of the product 
of eigenvalues. To find the integral over the amplitudes of the Fourier harmonics one has to 
express the amplitudes of the solutions chosen A k ,B k in terms of the amplitudes a k , b k . The 
relation between the coefficients A k ,B k and a k , bk can be found from the matching condition 
(first relation in 23) on the interface 



oo 



1—7= K cos(n9) + b n sin(n0)] = A g (w) _ _ [ A k^g k (w) + B k lmg k (w)} 

V 7r n=l T - Rs k=l 



r=R a 

(63) 

One can readily notice that the function g contains a large constant term, proportional to 
X, which totally dominates the matching condition for the ao mode at r = R s , so that 

-A ^[l + O(R/X)}. (64) 



v 7 ^ u 2(3 

For this reason the effect of the mixing between a and higher modes is suppressed by inverse 
powers of X and can be ignored in the limit of a long string. For this reason in considering 
the mixing of the modes in the following calculation we keep only k ^ 0. Furthermore the 
linear in x terms in the functions g and g\ are suppressed at r = R s by the factor R s /X 
and can also be disregarded. 

In what follows we consider the expansion of the functions g k at r = R s in powers of 
R s /{3, which expansion, as will be seen later, converges at R s < (3/2. Using 

oo 

(l + x)- k ^l + Y(-^ P C p p+k -^ (65) 
P =i 

where Cp +k _ 1 are the binomial coefficients, 

_ (p + fc-1)! 

Cp+k -' ~ p\(k-l)\ {bb) 
and also a definition of the Riemann ^-function 

oo 

C(A;) = £n-\ (67) 

we find for k ^ 

R k 00 / w \ p 
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with 



dpk 



-i) k + (-iy 



fj C( P + k)c; +k _ v 



(69) 



We have omitted in the expression (68) a constant term, which describes the mixing with 
the a mode as well as a term explicitly proportional to R s /X. Using the expansion (68) for 
gk{w) and considering the real part of the functions gk(w), one can find the coefficients a; in 
terms of Ak 

I p2tt 00 

ai = - V A k / Reg k cos(W)d6 = A x + V A k d lk , 
or in the matrix form 

A = {l + D)~ 1 a , 

where matrix D has elements dik- Similarly for the imaginary part one gets 



(70) 
(71) 



1 f 2w 

bi = - V fifc / lmg k cos{W)dO = B t - V B k d lk 
ii ^ Jo V 



and 



(72) 



(73) 



As usual, the contributions from the cos and sin modes are independent, and we consider 
the contribution to the boundary term from the even (cos) harmonics first 



R s 



de 



2tt 



de 



ai cos(w) 



L I 



E 

j k 



R \ k 00 
k[—j cos(kd) - pdpk 
" / P =i 



E^ 

k 



ka k - Y P a pd P k 



E a p i kS pk - pdpk] A k . 

k,p 



Introducing the matrix 



iV = diag(l,2,...,n,...) , 
one can rewrite the expression (74) in the matrix form 

Jo 



(R s g {R) d r g {R) ) \ r d6 = aN(l - D)A 



A substitution in this expression of the solution for A in terms of a (71) leads to 



2tt 



(R s gWd r gW 



d6 = aN(l-D)(l + D)- 1 a 



(74) 



(75) 



(76) 



(77) 
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Clearly, for the odd (sin) modes one gets the same expression with the replacement D — > —D. 
Collecting all the terms together one can write the result for the boundary partition functions 
(27) and (30) as 



y 12(boundary) 



Det 



(ei - e 2 )N 2 + e 2 N + e 1 N(l - D) 



1 I' 1 / 2 



1 + D 



{D - —D} 



%1 (boundary) — Det 



EiN + £iiV(l — D)- 



1 



'1/2 



{D -> -£>}, 



(78) 



l + £>. 

where {D — > — -D} means that one should take the expression and make the replacement 
D -> -£>. 

The zero temperature limit for the probability rate 7 S formally corresponds to setting 
D — > 0. Thus one can use the result for the zero temperature decay rate in Eq.(57), and 
concentrate on a calculation of the thermal catalysis factor K defined as 



dY 



£a Is 



(79) 



In a d- dimensional theory, i.e. with d — 2 transverse dimensions, the catalysis factor can be 
written as K, s = G d ~ 2 , where G is the factor per each transverse direction given by 



G 



7 crD=0 

^ / 12(boundary) l(boundary) 

5rD=0 ~ 

^ / 12(boundary) -^1 (boundary) 



(80) 



According to Eq.(78) it is a matter of simple algebra to express the factor G in terms of the 
matrix D: 

'N-l 



G = Det 



1 



N + b 



D 



'1/2 



(81) 



with b defined in Eq.(42). 



2.2.1 Analysis of the general formula 

In this section we consider in some detail the temperature effect in the string transition rate 
described by our general formula (81) in the situation where the inverse temperature is larger 
then the diameter of the classical configuration (bounce): (3 > 2R S . We first notice that due 
to the presence of the factor (N — 1) the first elements from the first row and the first column 
of the matrix D, d lk and d p \, enter the expression [(N — 1) D] 2 with zero coefficients, so that 
the final result (81) does not depend on them. Furthermore, one can see from Eq.(69) that 
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the matrix element d p k is not equal to zero only if the indices p and k have the same parity. 
Hence, there is no mixing between the amplitudes with even (02;, 62/) and odd (0,21+1, ^22+1) 
indices. Therefore the determinant in the (81) can be written as a product of determinants 
corresponding to the even and odd amplitudes 



Det 



1- 

\N + b , 



Det 



(2N-1 N 

1 - U 

\2N + b , 



Det 



1 - 



2 N 



2 JV + 1 + 6 



V 



(82) 



where the matrix elements of the matrices U and V are 



U n k — d 



Vnk — di 



pk — «2p2fc , Vpk — «2p+l,2fe+l , 



(83) 



and the indices p and k take values 1,2, ... . 

For practical calculations it is also convenient to write the expressions for the elements 
of the matrices entering in Eq.(82) in terms of their indices: 



XN + b U ) pk ~ 2v + b(2v)\(2k-lV a2p + 2k) [0) 



2N 



2N+l+b 



V 



pk 



2p + b (2p)\(2k-l)\ 

4p (2 P + 2k + iy (# 
2p+l + b(2p+l)\(2k)\ a2p + 2k + 2} {(3) 



2p+2fc+2 



,(84) 



with p, k — 1, 2, 

In order to find the first thermal correction at low temperature one can expand the 
matrices in power series using well known formula for the determinant 



Det(l — M) — exp [Tr In (1 - M)\ = exp 
In our case 



1=1 1 



1 — TrM + 0(M 2 



(85) 



( 36C 2 (4) (RX + 



2N - 1 
2N + b 



U 



600 C(6) 2 [RsY 2 120C(4)C(6) [R S \ W 



(2 + b) 2 \(3 (2 + b) (4 + 6) V P 



180C(4)C(6) (Rs 
(2 + b) (4 + b) {{3 



10 



(2 + by 

600 C(6) 2 



R.* 



12 



2iV 



2JV + 1 + 6 



2 / 1600 C(6) 2 ffl, 
I = 

V 



12 



(3 + 6) 2 ^ 



(2 + 6) (4 + 6) V/9 

\ 
• / 



(86) 
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Figure 8: The thermal catalysis factor per each transverse dimension vs RT for 6 = 1 (open 
circles) and for b = 10 (solid circles) 

Therefore the first correction to the zero temperature value of the rate is proportional to 
R^/P 8 and is given by Eq.(5). 

The series for the function G(T) diverges when (3 = 2R S . The corrections for the tem- 
perature close to (2R s )~ l can be found numerically. The proximity to this point defines the 
number of terms which should be taken into account in the series for G(T). The plot for the 
function G(T) vs the parameter R s /(3 = R S T calculated numerically with 50 first rows and 
columns retained in each of the matrices U and V is shown in Fig. 8. 

2.3 Decay of a string induced by two-particle collisions 

As we have seen in the previous sections, the only difference between the calculation of 
the decay at finite temperature T and at T = arises at the level of calculating the pre- 
exponential factor due to the functional determinant of the quadratic part of the action and 
amounts to the standard treatment of the boundary conditions in (Euclidean) time for the 
fluctuations: zero boundary conditions at large time for considering zero temperature and 
periodic boundary conditions with period (3 — 1/T at finite temperature. 

The formula (11) is in fact the unitarity relation [19] between the decay rate and the 
imaginary part of the amplitude of the transition amplitude from the false vacuum to the 
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false vacuum ((vacuum l) ou t ((vacuum l)i n )- Similarly, one can treat the probability of the 
string breakup by an excited state in terms of the bounce contribution to the imaginary part 
of its forward scattering amplitude. Proceeding to discussion of the string decay induced 
by the Goldstone bosons, we readily notice that a state with just one Goldstone boson 
cannot induce the destruction of the string. Indeed the total probability of such induced 
process is Lorentz invariant and thus can depend only on the (Lorentz) square of the particle 
momentum k 2 . The Goldstone bosons are massless, so that for them k 2 = and is fixed. If a 
single massless boson produced an effect on the decay, this effect would thus be not depend 
on the energy uj of the boson. In the limit uj — > the Goldstone boson is indistinguishable 
from the vacuum. (In other words the limit uj — > corresponds to an overall shift of the 
string in transverse direction.) Thus the decay rate of a single - boson state is the same as 
that of the vacuum, and the presence of a single boson with any energy produces no effect. 

The simplest excited state, contributing to the string destruction, is that with two parti- 
cles. The probability W 2 of creation of the critical gap in a collision of two particles with the 
two-momenta fa and k 2 can depend only on the Lorentz invariant s — (fa + k 2 ) 2 . Obviously, 
for two particles colliding on a string s = A.lo\L0 2 (and s = for two particles moving in 
the same direction, i.e. non-colliding). Using the unitarity relation, this probability can 
in principle be found in terms of the imaginary part of the forward scattering amplitude 
A(fa, k 2 ; fa, k 2 ): 

W 2 = c l ™ A (h,k 2 ;fa,k 2 ) 

UJ X UJ 2 

where the factor uj\uj 2 is the usual flux factor, and the constant C does not depend on 
either of the energies, and is determined by specific convention on the normalization of the 
amplitude. 

The dynamics of the Goldstone bosons on the string, including their scattering, can be 
considered in terms of the transverse shift Zi(x) treated as a two-dimensional field described 
by the Nambu-Goto action (8) as 

S = e 1 J A yjl + (d a Zi) 2 d 2 x + £2 J a yjl + (d a Zi) 2 d 2 x + /1 jf yjl + (dzi/dl) 2 dl , (89) 

where dl is the element of the length of the interface P between the phases. 

Clearly, at low energy of the Goldstone bosons one can make use of the expansion in 
Eq.(89) in powers of (dz) which generates the expansion of the scattering amplitudes in the 
momenta of the particles with each one entering the amplitude with (at least) one power of 
its energy uj, as is mandatory for Goldstone bosons. For the scattering in the metastable 
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state this generates expansion in powers of w/^/eT, so that in the lowest order in this ratio, 
that we are discussing here, it is sufficient to retain only the quadratic in (dz) terms in the 
action (89). It should be noted that in spite of retaining only the quadratic terms, the multi- 
boson scattering amplitudes do not vanish, since the necessary non-linearity arises from the 
bounce configuration. In other words, the bosons scatter 'through the bounce'. The energy 
expansion for these amplitudes is determined by the bounce scale £ c , so that the parameter 
of such expansion is (u£ c ) which we do not assume to be small. Clearly, the condition for 
applicability of the present approach where the terms of order ioj^fel are dropped, while 
those with the parameter oj£ c are retained is that 

^ (90) 



E\ — 62 S\ 

which is always true if the semiclassical tunneling can be applied at all to the string decay. 

We shall now show that in the on-shell scattering through the bounce each external leg 
enters with at least two powers of its energy. Let us start, for the simplicity of illustration, 
with the binary scattering. The general two — > two scattering amplitude A(k±, k 2 ; k 3 , k 4 ) 
can be related in the standard application of the reduction formula to the connected 4- 
point Green's function (vacuum 1|T \z (x i) z(x 2) z{xz)z{x 4)} (vacuum 1), which in turn is an 
analytical continuation of the Euclidean connected correlator 

(vacuum l\z(x 1 )z(x 2 )z(x 3 )z(x 4 ) (vacuum 1) = Z^ 1 — — — - | =0 . (91) 

5j{x 1 )5j{x2)5j(x 3 )5j(x 4 ) J 

The latter expression for the correlator assumes the conventional procedure of introducing 
in the action the source term / j(x) z(x) d 2 x for the Goldstone variable z(x) and Z b [j] is the 
path integral around the bounce in the presence of the source. 

The low-energy limit of the on-shell scattering amplitude is determined by the correlator 
at widely separated points Xi,...,x 4 . Weak 5-function sources 'prop' the string in the 
transverse direction at those points as shown in Fig. 9, and generally distort the bounce 
located between those points. The correlator (91) is determined by the distortion of the 
bounce by all four sources, so that at large separation between the points the bounce is 
located far (in the scale of its size £ c ) from the sources, where the overall distortion of the 
world sheet for the string is small and is slowly varying on the scale £ c . One can therefore 
expand the background field z s generated by the sources at the bounce location in the Taylor 
series around an arbitrarily chosen inside the bounce point x — 0. Clearly the first term 
of this expansion z s (0) corresponds to an overall transverse shift of the string and does 
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Figure 9: Bounce configuration distorted by the sources 



not alter the bounce shape and the action. Moreover the linear term in this expansion, 
proportional to the gradient dz s (0) does not change the bounce action either. Indeed this 
term corresponds to a linear incline of the string in the transverse coordinates, and can be 
eliminated by an overall rotation of the string in the d dimensional space. In other words 
the absence of the linear in the gradient of z s term in the action is a direct consequence of 
the d- dimensional Lorentz invariance of the string. We thus arrive at the conclusion that the 
expansion for the distortion of the bounce starts from the second order in the derivatives of 
z s (the curvature of the background world sheet), which for a connected correlator implies 
that in the expansion in the energy each external source enters with at least two powers of 
the energy. This conclusion clearly applies to the on-shell amplitudes with arbitrary number 
of external legs, since the generalization to multiparticle scattering is straightforward. 

In particular, the binary scattering amplitude A(k\, k 2 ', k 3 , k±) at low energy scales as the 
eighth power of the energy, and the related to it (by the Eq.(88)) probability of the string 
destruction by collision of two particles is proportional to the sixth power of the energy scale, 
or equivalently, to s 3 at small s. Applying in the same manner the unitarity condition to 
the forward scattering amplitude of a general n-particle state, one readily concludes that the 
corresponding probability W n of the induced string breakup scales with the overall energy 
scale uj as W n oc u 3n . 

2.3.1 Thermal decay and the string destruction by particle collisions 

The described procedure for calculating the bounce-induced scattering amplitudes in terms 
of the Euclidean correlators runs into the technical difficulty of calculating the bounce dis- 
tortion by the background created by the sources. Furthermore, this procedure obviously 
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involves a great redundancy, if the final purpose is to calculate the probabilities W n , i.e. 
only the absorptive parts of the forward scattering on-shell amplitudes are the quantities 
of interest. We find that in fact one can directly calculate the probabilities W n by an ap- 
propriate interpretation of the more readily calculable thermal decay rate in terms of the 
collision-induced probabilities. Such an approach is technically more tractable due to the 
fact that in the thermal calculation the bounce is not distorted, i.e. it is still a flat disk, and 
only the boundary conditions for the modes of fluctuations are modified. We first illustrate 
this approach by using the first nontrivial term of the low temperature expansion in Eq.(5) 
for calculation of the low energy limit of the binary probability W 2 . 

Indeed, the temperature dependent factors in the catalysis factor JC S arise from the con- 
tribution to the critical gap nucleation of the boson collision, weighed with the thermal 
number density distribution for the massless Goldstone bosons 

where k is the spatial momentum, \k\ = u, running from — oo to +00. Given the low energy 
behavior W n oc u 3n found in the previous section, one readily concludes that the contribution 
of n-particle string destruction starts with the term T 4n in the low T expansion for K. Thus 
the first term written in Eq.(5) can arise only from n = 2 i.e. from the binary production of 
the critical gap. Writing the expansion in s of the probability W 2 as W 2 = c 3 s 3 + . . ., one 
determines the coefficient C3 of the s 3 term by comparing the result in Eq.(5) with the one 
calculated in terms of the two-particle collision rate using the number density distribution 
(92): 



7T 8 f e 1 -e 2 \ 2 IIJT 
3ei-e 2 / V 2 



3 1R^6 



(d~2)c 3 r f 1 f T — 2 = {d-2)c^T\ (93) 
v ' 6 Jo e^/ T - 1 Jo e^' T - 1 4vr 2 v ; 3 225 ' K ' 

where the relation s = 4uiU2 is used and the integration is done only for the bosons with 
opposite signs of their momenta, and avoiding the double-counting for the identical particles. 
The factor (d — 2) counting the number of the transverse dimension, corresponds to the 
summation over the polarizations of the Goldstone bosons. The expression for the coefficient 
C3 following from Eq.(93) thus determines the first term in the expansion for W 2 



W 2 = ls R 



it 2 ( ei - e 2 x 2 



32 \3ei — e 2 



(94) 
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(One can notice that at E2 = the s 3 term coincides with the first term of expansion of the 
expression in Eq.(6).) 



2.3.2 Thermal bath with a chemical potential 

The discussed procedure for extracting the coefficients of the energy expansion for the proba- 
bility of collision- induced string decay is obviously limited to only the first term in W 2 . In the 
higher terms in the temperature expansion for K, s the contribution of the energy expansion 
for W n with different n generally gets entangled. This happens because the temperature is 
the only parameter and terms originating from different n can contribute in the same power 
in T. In order to disentangle the terms of higher order in energy in W n with low n from 
similar terms originating from higher n we introduce a negative chemical potential v for 
the Goldstone bosons. Generally such procedure would be impossible, since the number of 
these bosons is not conserved. However in our application this procedure is fully legitimate. 
Indeed the thermal state of the string that we study here is that of collisionless bosons, in 
which their number is conserved. The string decay, resulting in a change in this number, is a 
very weak process that we consider only in the first order, which justifies averaging the rate 
of this process over the unperturbed state with conserved number of particles. At negative 
v the number density distribution of the bosons (92) is replaced by 

dn{k) = -^ ^, (95) 



e t 



I 2tt ' 



and by tuning the parameter \v\jT one can readily resolve between the contribution of 
n-particle processes with different n. 

The introduction of the chemical potential requires us to modify our previous thermal 
calculation (see Sec. 2.2). The modification of this calculation for a thermal state with a 
negative chemical potential, where the number density of the bosons be given by Eq.(95), is 
achieved by introducing a 'damping factor' in the periodic sums for the outer functions in 
Eq.(61): 



X Dfe Dfe 

w K — npp (w + np)" 



exp(—n\v\/3) . (96) 



One can readily see that the only net result of the v dependent factor in the calculation is a 
modification of the matrix coefficients d p k amounting to the replacement of the factors ((q) 
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by the standard polylogarithm function, 

oo 

Li q (x) = 



n=l 



X'" 



C(q) — > Lig(e ^^ T ). In other words, the catalysis factor for a thermal state with a negative 
chemical potential is given by the expression 

K 8 (v, T) = Det [l - V 2 (u, T)J Hd ' 2) ' 2 

with the elements of the matrix T>{y,T) having the form 

,k+p+2 P (p + k + 1) 



(97) 



P(v,T)} pk = -[(-iy + (-iy\ {r s t) 



p + b (p+l)\k\ 



(98) 



The dependence on two 'tunable' parameters v and T in Eq.(97) makes it possible to 
disentangle the contribution of processes with different number of particles from the energy 
behavior in each of these processes. Such a separation becomes straightforward if one notices 
that each factor with the polylogarithm Li arises from the integration over the distribution 
(95): 

(99) 



f 

Jo 



00 H 



e t —l 



du = q\ T q+1 Li q+1 (e 



W\/t 



One therefore concludes that the number of the l Li factors' in each term of the expansion of 
the catalysis factor in Eq. (97) directly gives the number of particles in the process, while the 
indices of these polylogarithmic factors give the power of the energy for each particle. Given 
that the matrix T>(u,T) is linear in the l Li factors', one can count the number of particles 
contributing to each term of the expansion for K,{y,T) s by counting instead the power of 
T>(u,T). The latter counting is simplified if one rewrites the equation (97) in the equivalent 
form, suitable for the expansion in powers of V(u,T): 



1 + 



Tr 



K, s (u, T) 



V(y,Tf 



+ 



exp 
d- 



d-2 



Tr In 



Tr 



V(u,T) 4 



+ 



8 



(100) 



{Tr [v(u,T) 2 ]} 2 + O (V 6 ) . 



The latter expression merits some observations. The first is that all the terms in the 
expansion in powers of V(u,T) are positive, which is certainly the necessary condition for 
the consistency of our interpretation of these terms as corresponding to the probability 
of the destruction of the string by n-particle collisions. The second is that the string is 
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destroyed only in collisions of even number of particles, since the expansion goes in the even 
powers of T>(u, T). Finally, the third observation is related to the dependence in Eq.(lOO) on 
the number of the transverse dimensions (d — 2). Namely, the quadratic in V(u,T) term is 
proportional to (d— 2). This corresponds to that in two-particle collisions only the Goldstone 
bosons with the same transverse polarization do destroy the string. On the contrary, the 
quartic in T>(u,T) term has one contribution proportional to (d — 2) and one proportional 
to (d — 2) 2 . The linear in (d — 2) part corresponds to all the bosons in the collision having 
the same polarization, while the quadratic in (d — 2) part is necessarily contributed by the 
collisions, where the colliding bosons have different polarization. 

2.3.3 Destruction of the string in two-particle collisions at arbitrary R s ^/s 

The expression (100) for the catalysis factor JC s (v,T) together with the formulas (98) and 
(99) reduce the calculation of the probability of the string breakup by a collision of an 
arbitrary (even) number n of particles to straightforward, although not necessarily short, 
algebraic manipulations. In this section we consider in full the most physically transparent 
case of two-particle collisions. The probability in this case is found from the term in Eq.(lOO) 
with the trace Tr[D(v,T) 2 ]. Using Eq.(98), this trace can be written as a double sum: 



Tr 



V(u,T) 2 ] = (101) 

(p + k + 1)! (R s T f +k+2 U p+k+2 (e-M/^)] 2 



4 §£ ( P + b + l)(k + b+l) (p-l)!(p+l)!(*-l)!(* + l)! 

One can readily recognize the expression in the straight braces here as the integral (99) with 
the power of the energy q given by (p + k + 1), and thus identify the coefficient of the same 
power of s = 4ujiU2 in the expansion of the probability W^(s) in s. In this way we find the 
following formula for W^(s) in terms of this expansion, 

oo oo n / r>2/A\p+ k + 1 

W 2 (s) = 8vr 2 7s ^£E {SRJA) 



tit^i (p + b+l)(k + b+l) (p-l)!(p+l)!(A;-l)!(A; + l)! 
= 8tt 2 7sj R s 2 [MV^Rs)f , (102) 

where the function $f>(x) expands in a single series as 

*>(*Hfg p+ L 1(y -i)f; + i)r « 
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It can be noted that the two-particle probability depends only on the odd powers of s. This 
in fact is consequence of the binary forward scattering amplitude being even in s, as required 
by the Bose symmetry. 

For integer values of the parameter 6, the function $b( x ) has a simple expression in terms 
of the modified Bessel functions I q (x) of the order q up to q = 6 + 2. This expression is 
especially simple for 6 = 1, i.e. for the case of the string decay into 'nothing': $i(x) = 
Is(x), so that one arrives at the formula (6). We also write here, for an illustration, the 
corresponding expressions for the next two integer values of 6: 

$ 2 (z) = - [h(x) + 6I 4 (x)} ; $ 3 (z) = \ [(48 + x 2 ) h{x) +xh{x)] . (104) 

Ob Ob 

3 Domain wall 

3.1 Spontaneous decay of a metastable domain wall 

In the following sections we generalize the previously obtained results for the decay of a 
metastable domain wall. We denote the tension of a domain wall by e and the tension of a 
string associated with the edge of the wall by a, so that there should not be any confusion 
with previously used notations for a string. For a complete decay of a domain wall the low 
energy effective action is similar to (8) and is given by Nambu-Goto action for two and tree 
dimensional objects 

S = aA + eV, (105) 

with V being the world volume of the wall, while A is the world area of interface. As before 

this action does not take into account the inner structure of the wall or the interface. A 

nontrivial classical solution (bounce) in this case is empty sphere with the radius 

2(7 , . 

R w = — , (106) 

surrounded by the metastable phase (Fig. 10). The appropriate quantity to be found is j w : 
the nucleation rate of critical holes per unit area (not length as it was for a string transition). 
Following the same steps as in calculating the decay rate of a string we find 2 

7. = ^expf-^-^l — , (107) 



2 The factor exp {— 16 3 ^ ^ corresponds to the decay rate of a metastable vacuum in 2 + 1 dimensions 
[12]. 
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Figure 10: Bounce configuration 



with boundary partition functions given by 



^(boundary) = f V£ eX P ~ J dft£A«£+ ^ J d{l£d r Z, 



r=Ru 



and 



Z\ (boundary) = J T>£ exp J 



D2 

w 1 dn^d r z t 



r=R 



'108) 



(109) 



with the function z c (r, 9,ip) satisfying the Laplace equation Az c = with the boundary 
conditions 

z c (R w ,0,<p) = £(0,<p), z c (r = L) = 0. (110) 

The operator A( 2 ) is the angular part of the Laplace operator in 3d (the Laplace operator 
on a sphere). Expanding the boundary function <p) in the series in spherical harmonics 



we find a complete set of functions 



111) 



Lm 



Z lm (r, 9, if) = Airn—^- Yl m (9, if) . 



R l+1 



(112) 
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Substituting these functions to the equations (108) and (109) and performing integration 
over the amplitudes A im yields 



'12(boundary) 



A/-n 



(2/+l)/2 



and 



f-J \al(l + l) + eR w (l + l)J 

» / 1 \ (21+1) /2 



-2^1 (boundary) — A/" J| 



^ Vei^ (Z + l). 



(113) 



(114) 



3.1.1 Regularization 



Each of the expressions (108) and (109) contains a divergent product (compare to Eqs.(34) 
and (35)). As previously, introducing Pauli-Villars regulators with masses M Q , we find the 
ratio of the boundary partition functions to be equal to 
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1=0 
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1=0 
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£=0 



/(/ + 1) + 2JM^2 +(/ + !) +1 



(/ + !)(/ + 2) 



/ + 1 



(2Z+l)C a /2 



X 



(21+1) C a /2 



M 2 a Rl+(l + l) +1 



1 + 



x 



MlRl 



MlRl+(} + \) 2 ) (l{l + l) + 2sjMlRl+(l + \) 2 + l^ 



(21+1) C a /2 

(H5) 



where we took into account that R w = 2cr/e. 

Now, each of the products in Eq.(41) is finite at a finite M a and can be calculated 
separately. Instead of calculating the product directly, we can use the relation 



(116) 



and calculate the sum. We start from the third product. The expression under the sign of 
product is of the form 1 +g(l), with g(l) close to for any /, since it behaves as M~ l . Hence 
we can leave only the first term in the expansion of the logarithm ln(l + x) = x + 0(x 2 ). 
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Thus, we need to find the sum 

oo 



1 



s» = ^ EE Ca (2m) mX 



M 2 i? 2 + [/ + 



+ Z(l + l) + 2^i22,+ (z + i) +1 



(117)" 

The sums associated with the three products are readily calculable with the help of the 
Euler-Maclaurin summation formula and the result for the regularized ratio has the form 



1Z W = exp 



M 2 R 2 In MR W + MR W In MR W + In MR V 



(118) 



where M n lnM = J2 a C a M™ In M a , for any n. The expression for 1Z W contains an essential 
dependence on the regulator mass parameter M. However, similarly to the previously consid- 
ered case all such dependence in the phase transition rate can be absorbed in renormalization 
of the parameter a in the leading semiclassical term. 

Using the same technique as employed in Section 2 one can find the renormalized param- 
eter a to be given by 

a R = a + 5a 



(119) 



with 
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(120) 



3.1.2 Results 



Collecting all terms together we find the rate of the process. It is clear that the result for 
each d — 3 transverse dimensions factorizes, thus we have the expression for the rate in the 
form 

;i2i) 



7« 



Tit 3 exp 



r7/3 



3e 2 



where a is the bare (non-renormalized) tension, and the regularized ratio 1Z is given by (118). 
Taking into account that each of the transverse dimensions contributes additively to 5a and 
the interface is a sphere with area 

A = 4nR 2 w , (122) 
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and expressing the bare er through the renormalized one: a = <jr — 5a, one readily finds that 
the dependence on the regulator mass M cancels in the transition rate, and one arrives at 
the formula given by Eq.(4). 

Thus we obtained the result similar to the decay of a string, where the effect of all extra 
transverse dimensions results only in the renormalization of parameter a associated with the 
interface. It should be mentioned, however, that the result for a string was, actually, obtained 
for a transition between two states of a string with different tensions. Here we considered 
decay of a wall (transition into nothing). For the calculation used it is impossible to preserve 
finite terms, but only proportional to some power of the regulator mass parameter M. 

Having calculated the probability rate for a decay of one- and two- dimensional objects, 
e.g. string and domain wall, it is tempting to assume that the situation is somewhat similar 
for the decay of an object of arbitrary dimensionality. But it is not true already for the decay 
of tree- and four- dimensional objects, where the dependence of the result on regulator mass 
is substantial even after renormalization of a parameter associated with an interface. That 
dependence demands introduction of new terms into initial action, which corresponds to 
nonrenormalizibility of the effective 'low-energy' theory. 

3.2 Thermally induced decay of a domain wall 

The thermal effects in decay of a metastable domain wall involve one important difference 
from those in the decay of a string. Namely, the previously considered expansion of the catal- 
ysis factor K s is valid up to the temperature T = 1/ i c = 1/(2R S ), at which point the whole 
calculation breaks down due to a change in the mechanism of the transition: the thermal 
fluctuations of the string start dominating over the quantum ones. Simultaneously, at this 
temperature the thermal factor in Eq.(81) develops a singularity, and the low-temperature 
expansion diverges. In the case of a domain wall decay at low temperature the thermal 
effects in the quantum tunneling result in a catalysis factor K w multiplying the semiclassical 
exponential factor in •y w : 



On the other hand, the static potential for a bubble, a round hole in the wall, depends on 
the radius R of the bubble as 



(123) 



V(R) = lixaR-ntR 2 , 



(124) 
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and this potential has a maximum at R m = a/e where its value is 



ix a 



(125) 



The probability of classical thermal fluctuations over the barrier V m is proportional to the 
activation factor exp(— V m /T), and this factor becomes larger than the exponential term 
in ^ w (Eq.(4)) starting from T = 3/(8R w ), at which temperature the classical fluctuations 
over the barrier start to dominate over the considered here quantum tunneling. In terms 
of the Euclidean space formulation of the problem the periodic replication with the period 
/3 of the spherical bounce (Fig. 11a) gives a larger action per period[15] than the cylindrical 
configuration in Fig. lib as soon as T > 3/(8R w ). For this reason our calculation of the 
thermal factor K, w makes physical sense only as long as the temperature is lower than this 
value. 
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(a) (b) 

Figure 11: Bounce configuration for the quantum tunneling at low temperature (a) and for 
the thermally activated transition (b). 



In calculating the preexponential factor arising in the periodically replicated bounce 
configuration, we encounter the same problem of a mismatch between the symmetry of the 
boundary conditions and the symmetry of the bounce. However one can cope with this 
problem in a way similar to the previously discussed approach to calculating the decay rate 
of the string at finite temperature. Namely, one can consider the solution to the Laplace 
equation as one produced by the sources placed at each period (see Fig. 12). Thus the 
solution Z[ m (the periodic analog of the functions in Eq.(112)) is given by 

R i+i R i+i 

Zlm = -gr*W0, <P) + E -TT Y lrn(0n, <f), (126) 
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with r n and 6 n reading as 



r n = \Jr 2 + (f3n) 2 — 2(3nr cos 6, 

. r, rsm.6 /,^„x 
sin# n = . (127) 



The solution (126) can be expanded in a power series 
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Figure 12: The illustration of the notation for the angles and distances used in Eq.(127). 



In order to find the constants CJJj one can notice that the expression in Eq.(126) has the 
same structure for each n. Therefore, for the purposes at hand we can consider only two 
adjacent periods (±/3). Then integrating the r.h.s. of the (126) with Y^(0, </?) we find the 
result to be 



qj? = 2(-l) l+l ' +m J^±A + (129) 

V 21 + 1 J {l + m)\ (/ - m)\ {V + m)\ (/' - m)\ 



It should be mentioned here that obviously 



E <T = 0. (130) 

m=—l 
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This relation can be understood even without actual calculation. Indeed, the expression 

/ J2 Y i*Mv)Yim(0i,<p) sine de dip (131) 

is proportional to the electric potential at point (3 created by the set of 2'-poles placed on a 
sphere with its center in the origin. Obviously, such a potential is equal to zero. 

Once we have found the coefficients (129) of the expansion (128) we can follow the same 
routine as used in the section (2.2). Namely, we can express any solution of the Laplace 
equation with periodic boundary conditions as a series 

z c = ^zi m Bi m . (132) 

Ira 

After that using the expansion of a boundary function (111) and the relation (110), we 
express coefficients Bi m through A[ m 



Aim — 



fe+fyj C%C(l + l' + l) 



Bi'm, (133) 



or 

A m = (l + D m )B m , (134) 
where D m is a matrix with elements 

[D m ]w=l^jj C%((1 + 1> + 1), (135) 
and A m (similarly B m ) is the row type object with elements 

[An]l = Aim- (136) 

It should be mentioned that since the solution z\ m exists only for I > m, the row A m has 
nonzero elements A\ m only for I > m. 

Substituting the solution (132) to the expressions for the partition functions (108) and 
(109) and integrating over the amplitudes A irn (with a help of relation (134)) we find the 
catalysis factor in the following form 

k -^ a A x +tMwt2) d 'V ' (137) 
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with matrix N defined in (75). Equivalently, one can rewrite the expression in the following 
form 

_ d-3 

K w = det (l + £>) 2 , (138) 
where V stands for the block diagonal matrix defined as 



V 



\ 



(139) 



with elements 



N(N-l) 



Dr 



(N+l)(N + 2)~ m ' (140) 
The first nontrivial term of expansion of the catalysis factor JC W at small temperatures is 
determined by the coefficients in Eq. (135) at / = /' = 2 and the result for this term is given 
in Eq.(7). 

One can notice that the low-temperature expansion for the factor K, w generated by 
Eq.(138) is well behaved at TR W = 3/8, beyond which temperature, as previously dis- 
cussed, the regime of the transition changes from quantum tunneling to a classical thermal 
activation. Moreover, the calculated thermal effect at this point is numerically extremely 
small: K, w — 1 ~ 10~ 3 . Thus the explicit form of higher terms in the temperature expansion 
is of a 'practical' interest only inasmuch as the thermal calculation is used as a preliminary 
step for finding the generating function for the rate of the wall destruction in collisions of 
Goldstone bosons. 



3.3 Destruction of a metastable domain wall in binary collisions 

Using the arguments from the section 2.3.2 it is possible to relate the catalysis factor K. w for 
a decay rate of the domain wall at finite temperature calculated in the previous section, to 
the effective length of a particle collision (the analog of the probability in 1 + 1 dimensional 
case) . 

As before we consider the distribution function for the Goldstone bosons in the following 
form (with negative chemical potential — \u\) 

1 d 2 k 



dn(k) = 
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(141) 
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Such an introduction of the chemical potential modifies the result for the catalysis factor 1C W 
in the way that the solutions (128) are modified as 

z lm = ^^) + EE(^^(^)^ (142) 

JDl+l _ T>1' l+l 

= ^T^(^' ^) + E + fc + ^^'"(^ ^) Li Hi' + i W )- (143) 

This modification can be accounted for by the substitution £(/ + I' + 1) — > Li2 + p + i(e~'^/ T ). 
Thus, the catalysis factor becomes (with d — 3 transverse dimensions) 



JC W = det 



d-A 
2 



1 + P(7» 2 , (144) 



where the matrix P(T, z/) has a block diagonal form with elements 

[Vm]w(T, v) = y + l^ly C w (RT) W ' +1 Li w+1 (e-l^). (145) 

In order to separate the contribution of binary collisions we expand the catalysis factor 
in powers of V to the second order: 

K w = 1 + ^TrV\T,u) = ^ETr^, (146) 



m 



where we have taken into account that the linear term vanishes: 

TrP(T,z/) = ETrP m = 0, (147) 

m 

due to the relation (130). Similarly to the case of string decay the physical reason for the 
absence of the first power is the Lorentz invariance of the system (there is no destruction 
of the wall by presence of one massless particle). It should be also noted that unlike in the 
decay of a string the destruction of a domain wall can occur in collisions of an odd number of 
particles, since the expansion (146) generally contains any larger then one power of matrix 
V. 

For an actual calculation of binary processes one can now follow the same steps as in 
the calculation of the section 2.3.2. The only difference from the (1+1) dimensional string 
geometry is that now the kinematical invariant s also depends on the angle 6 between the 
two particles' momenta, 

s(k 1: k 2 ) =2cj 1 cj 2 (1 + cos^) . (148) 
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Therefore an analog of the integral (99) has an angular part, and one should make use of 
the relation 

d 2 h d 2 k 2 4™" 1 



/ n(h) n{k 2 )s-\K £ 2 )^P = ^ (n!) 2 T 2 ^ Li 2 l+1 (e^) B(n - 1/2, 1/2), 

(149) 

where B(a, b) is the Euler beta function 



(160) 

The appropriate quantity describing the probability of a binary process on a (2+1) dimen- 
sional domain wall is the effective length A, which is an analog of the cross section in (3+1) 
dimensions and of the dimensionless probability W 2 on a (1+1) dimensional string. 

Using the expression (149) and the quoted in Eq.(7) low temperature behavior of the 
thermal catalysis factor one readily finds the low energy limit for the effective length of 
destruction of the wall in a collision of two Goldstone bosons: 

X = ^7r 2 lw s 3 R 1 w ° + ... . (151) 

The general formula for the effective length at arbitrary values of y/s R w is found using 
Eq.(146) and the expression (145) for the elements of the matrix T>. In this way we arrive 
at the following expression for the effective length of a two particle collision in the form of a 
triple sum 

A = - 1671-3 (— ~ 3 ) R ^ lw J2 viv-^pip- 1 ) 



t%,^ m |(? + l)(? + 2)(p + l)(p + 2) 



/ D fZ\ 2p+2q I 

x £WM B -if p + q _ 1/2,1/2)- — -f- — r7,(152) 

where j w is a decay rate per unit area of the wall at zero temperature. 

The behavior at large energy, R w y/s 3> 1, can be found using saddle point approximation 
to estimate the sum over p and q, which gives 

A~expf — +4^-1, (153) 

which matches the semiclassical expression for the tunneling exponent at energy y/s [22]. 
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4 Summary and conclusions 



In this paper we considered the decays of metastable topological configurations such as string 
and domain wall. For each process we found the rate at zero temperature and calculated 
catalysis factor at finite temperature. Using the relation between catalysis factor and proba- 
bility (effective) length of a collision of the Goldstone bosons we found probability of a string 
(domain wall) decay in a collision of two particles. The main conclusions that can be drawn 
from the presented here study can be formulated as follows 

i. The effect on the preexponential factor in the tunneling decay rate arising from the 
Goldstone modes living on strings and walls is expressed in terms of boundary terms in 
the path integral around the bounce configuration. This effect is fully calculable within 
the low-energy effective action of the Nambu-Goto type, and essentially reduces to a 
renormalization of the tension of the interface between the two phases of the topological 
defect. We have found that this property is unique for the strings and walls and is 
generally lost in a description of decay of topological objects with more dimensions, 
where one would have to go beyond the effective Nambu-Goto description. 

ii. The thermal effects in the tunneling rate are also fully calculable and are described by 
an expansion in powers of the temperature. The power-like, rather than exponential, 
behavior of the thermal terms is due to the presence of the massless Goldstone bosons. 

iii. The expansion for the thermal catalysis factor converges as long as the temperature 
is less than the inverse diameter of the bounce, so that in the Euclidean calculation 
of the bounce configuration to (the imaginary part of) the free energy there is no 
obstruction to periodically replicating the bounce with the Matsubara period (3=1 /T. 
For the decay of strings this range of temperatures corresponds to the dominance of the 
(thermally enhanced) quantum tunneling over the classical thermal effects, and at the 
critical temperature l/£ c the regime of the transition changes and the classical thermal 
effects take over. For the domain wall decay the classical thermal activation starts 
dominating at a lower temperature, where the thermal enhancement of the quantum 
tunneling is still very small. 

iv. The thermal enhancement of the tunneling rate can be viewed as and additional con- 
tribution to the decay probability arising from collisions of the Goldstone bosons which 
are present in the thermal state. Such an interpretation in fact allows to determine the 
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rate of destruction of the considered topological objects in collisions of the Goldstone 
bosons. In order to unambiguously identify in the thermal expansion the terms corre- 
sponding to the contribution of collisions between a given number of the bosons, we 
have modified the considered statistical ensemble by introducing a (fictitious) negative 
chemical potential v for the Goldstone bosons. The calculated dependence on v and 
T of the preexponential factor in the decay rate thus serves as a generating function 
for the rates of the decay induced by collisions of particles, including the dependence 
of those rates on the energies in the collisions. 

v. We find that a destruction of a string takes place only in collisions of even number 
of the Goldstone bosons, while a metastable wall can be destroyed in collisions of any 
number of the bosons, n > 2. 

vi. We have calculated in a closed form the energy behavior of the probability of destruc- 
tion of a string in a collision of two Goldstone bosons, which is valid at an arbitrary 
relation between the energy ^/s and the infrared scale in the problem R s . The found 
expression takes an especially simple form of Eq.(6) for the decay of a string into 
'nothing' 

vii. At large values of R -Js the collision-induced decay rate contains an exponential factor 
exp(2_R v /i), which has been previously argued[17, 22] within the leading semiclassical 
approach and corresponds to a full transfer of the energy from the colliding particles 
to the tunneling degrees of freedom. 

One can readily notice that our calculation of the collision-induced rate through a thermal 
ensemble with an artificially introduced chemical potential is quite indirect. In a sense, this 
approach is reminiscent of the treatment in Ref. [23], where finite volume effects in the 
energy of a two-particle state are related to the binary scattering amplitude. The obvious 
difference is that we calculate the free energy of a statistical ensemble rather than of a 
particular quantum state. In relation to this part of our calculation it would certainly 
be illuminating to have a more direct, and possibly simpler method for calculating the 
probability of creating semiclassical objects, such as the bubbles of the stable phase, in 
collisions of particle. However, lacking such a method at present, we have to resort to the 
indirect calculation described in the present paper. 
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